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Realization  of  Metamaterial-Based  Devices: 

Mathematical  theory  and  Physical  Demonstration 

Hossein  Mosallaei 
ECE  Department 
Northeastern  University 
Email:  hosseinm(3iece.neu.edu 

This  project  focuses  on  dipole  mode  analysis  and  dispersion  diagram  characteristics  of  metamaterials 
constructed  from  array  of  dielectric  and  magnetic  spheres  in  a  unit-cell.  A  comprehensive  study  is 
provided  and  different  metamaterials  are  established.  It  is  demonstrated  combining  electric  and 
magnetic  dipole  modes  (different  spheres  particles  in  a  unit-cell)  can  successfully  engineer  the 
dispersion  diagrams.  Effective  metamaterials  parameters  are  obtained.  Bandwidth  enhancement  by 
increasing  the  coupling  between  the  particles  is  also  realized.  The  performance  analysis  of  other 
shape  structures  (such  as  disks,  easier  for  fabrication)  is  also  explored.  At  the  end,  hemisphere 
magneto-dielectric,  ENG,  and  MNG  small  antennas  are  developed  and  their  physics  are  highlighted. 
Our  theoretical  model  is  successfully  validated  using  a  full-wave  numerical  approach  through 
successful  modeling  of  complex  metamaterial  structures.  It  must  be  mentioned  that  in  this  report  a 
detailed-summary  is  presented,  where  the  full  and  comprehensive  discussions  can  be  found  in 
published  papers  and  dissertation  made  in-part  out  of  this  work  (see  references).  One  student  has 
been  fully  working  on  this  project  (however  some  of  my  other  students  have  also  helped  with 
discussions  and  some  simulations) 


1.  Array  of  Spheres 

Modeling  of  array  of  dielectric  spheres  (see  Fig.  1),  with  arbitrary  geometries  and  locations,  are 
explored.  The  incoming  and  scattered  fields  for  each  sphere  is  expressed  in  terms  of  Mie-series  as 
below  [1] 
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And  then  using  the  additional  theorems  and  applying  the  boundary  conditions  the  required  unknown 
parameters  are  derived: 
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Where  the  summation  over  q  is  the  summation  on  all  particles  except  the  reference  sphere,  is 
the  center  distance  between  p""  sphere  and  q"'  one.  and  are  the  translational  addition 
theorem  coefficients.  Equations  (3)  are  a  coupled  set  of  complex  linear  algebraic  equations,  and 


should  be  solved  simultaneously  to  yield  the  unknown  scattering  coefficients.  In  addition,  the  infinite 
series  must  be  truncated  to  a  finite  number  n  =  v  =  N  and  m  =  /u  =  M.  To  obtain  the  scattering 
coefficient  first  we  have  to  determine  the  vector  translation  additional  coefficients.  Based  on  the 
Cruzan's  theorem,  the  translation  addition  coefficients  are  [2]: 
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In  the  above  equations  represents  the  Hankel  function  of  the  first  kind,  P™  is  the  associated 
Legendre  function  of  the  first  kind.  The  a(-m,n,^,v,p)  and  a{-m,n,/j,v,p,p-l)  are  the  Gaunt 
coefficients  [3].  The  Gaunt  coefficient  are  closely  related  to  Wigner  3jm  symbol  and  defined  by: 
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The  other  coefficients  in  (3)  are  defined  as; 
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Note  that,  the  integer  p  in  the  summation  of  (3a)  takes  the  values  of  |r)-v|,|r)-v^|+2,..., n  +  v' .  While 
this  integer  (p )  in  (3b)  varies  from  |  n-v  |+1,  |r)-v  |  +3,...,  a7  + v  .  The  Wigner  3jm  symbol  is  defined  as: 
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Once  we  calculated  all  the  addition  theorems,  the  scattering  coefficients  can  be  calculated  using  (3). 
Now  that  we  have  the  scattering  coefficients,  the  total  scattering  field  of  the  entire  cluster  can  be 
obtained.  In  the  far-field  region  the  relation  between  the  electric  vector  components  of  the  incident 
field  and  those  of  the  scattered  filed  is  given  by; 
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where  k  is  the  wave  number  and  r  is  the  usual  spherical  coordinate  for  a  point  in  space.  The  electric 
vector  components  (f  ||,,f  1,)  and  of  the  incident  and  scattered  fields  are  parallel  and 

perpendicular  respectively  to  the  scattering  plane  defined  by  the  direction  of  the  incident  wave 
vector  and  the  scattering  direction.  To  check  the  accuracy  of  the  scattering  coefficient,  we  compare 
the  theoretical  results  of  components  of  polarized  scattered  intensity  /‘n  and  /jj  with  the 
experimental  results  in  [4,  5].  /„  and  are  functions  of  the  scattering  angle  0  alone, when  referred 
to  the  scattering  plane  only  (^  =  0)  and  can  be  calculated  from  the  amplitude  scattering  matrix 
through  the  equations 

/„=|sr(0)P  (9a) 
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where  the  superscript  x  or  y  indicates  that  the  incident  plane  wave  is  x  or  y  polarized  (note  that  the 
incident  plane  wave  propagates  along  the  z  axis).  For  each  sphere,  and  are  defined  as  follows: 
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The  scattering  parameters  for  two  types  of  structures  are  obtained  in  Figs.  2  and  3.  The  results  are 
compared  very  successfully  with  experimental  demonstrations  in  [4,  5],  validating  the  accuracy  of  our 
technique.  Note  that,  for  all  the  cases,  the  axis  of  symmetry  of  each  sphere  chain  is  perpendicular  to 
either  the  scattering  plane  x-z  plane  or  the  incident  wave  vector  along  the  z  axis.  In  other  words,  it  is 
always  parallel  to  the  y  axis. 


z 


X 

Fig.  1:  Configuration  of  an  aggregate  of  spheres. 


Fig.  2:  Comparison  of  theoretical  calculations  with  laboratory  scattering  measurements  from  [5]  for  angular 
distributions  of  polarization  components  of  scattered  intensity  by  two  identical  spheres  with  inter-sphere 
distance  of  10.76.  Each  sphere  has  a  size  parameter  of  ka=4.346  and  a  refractive  index  of  er=2.6568  +  0.0326i, 
(a)  log  (ill),  (b)log(i22). 


(a)  (b) 

Fig.  3:  Comparison  of  theoretical  calculations  with  laboratory  scattering  measurements  from  [5]  for  angular 
distributions  of  polarization  components  of  scattered  intensity  by  two  identical  spheres  with  inter-sphere 
distance  of  9.94.  Each  sphere  has  a  size  parameter  of  ka=4.346  and  a  refractive  index  of  er=2.6568+0.0326i, 
Hr=l  (a)  log  (ill),  (b)  log(i22). 


2.  Two-Set  Dielectric  Spheres  Unit-Cell 

The  objective  in  this  section  is  to  investigate  metamaterials  constructed  from  arrays  of  two-set  of 
spheres  in  a  unit  cell,  as  illustrated  in  Fig.  4.  Considering  one  is  exploring  a  metamaterial  realization, 
the  spheres  sizes  must  be  subwavelength  and  thus  one  can  use  the  equations  derived  in  section  1 
and  simplify  them  by  considering  only  the  electric  and  magnetic  dipole  modes.  This  will  extensively 
reduce  the  complexity  of  the  modeling.  Figs.  5(a)  and  5(b)  demonstrate  the  scattering  performance  of 
a  single  sphere  around  the  frequency  band  that  only  the  electric  and  magnetic  modes  exist.  As  one 
can  see  the  scattering  coefficients  goes  to  resonance  but  at  different  frequencies.  To  achieve  a 
metamaterial  with  tailored  electric  and  magnetic  parameters,  one  needs  to  bring  these  two 
resonances  around  the  same  spectrum.  This  can  simply  be  achieved  by  making  a  metamaterial 
constructed  from  two  spheres  in  a  unit  cell  (the  same  materials  and  different  sizes  or  the  same  sizes 
and  different  materials). 


The  structure  is  a  complex  configuration  and  one  needs  to  comprehensively  model  the  metamaterial. 
A  reciprocity  technique  is  applied  that  the  fields  on  a  reference  sphere  is  equal  to  the  scattered  fields 
from  all  other  spheres  at  that  element.  Using  the  periodic  boundary  condition  and  the  fact  that  a 
travelling  wave  with  proper  phase  shift  exists  through  the  array,  the  below  equation  for  achieving  the 
dispersion  diagram  for  metamaterial  (kd-pd)  is  achieved  [1,  6]: 
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5”!  and  S[  for  r€{l,2}  are  the  normalized  Mie  electric  and  magnetic  scattering  coefficients. 
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g^(a)  and  ^j(o)  are  the  first  Mie  electric  and  magnetic  scattering  coefficient: 
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Here  is  the  wavenumber  inside  the  p"’  sphere  while  and  are  the  spherical  Bessel  and 
spherical  Hankel  function  of  the  second  kind.  The  summations  defined  in  (11)  are  given  by  [6]: 
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It  is  worth  noticing  that,  the  summations  in  (14)  are  the  interactions  between  the  scattered  electric 
(magnetic)  field  from  other  spheres  and  the  electric  (magnetic)  field  of  reference  spheres  while  the 
terms  in  (15)  are  the  coupling  between  the  scattered  electric  (magnetic)  field  from  other  spheres  and 
the  magnetic  (electric)  field  of  reference  spheres.  Also,  one  can  divide  each  of  the  series  in  (14),  and 


(15)  into  two  parts,  one  the  contributions  from  the  self  plane  spheres  (/,  =0)  and  other  the 
contribution  form  all  others.  For  instance,  can  be  expressed  as; 

(16) 

Where  is  the  summation  over  the  scattered  field  from  spheres  in  the  self-plane  (the  plane  of  the 
reference  sphere)  and  is  the  summation  of  the  scattered  field  over  all  other  planes  except  l,=0. 
These  expressions  help  us  better  understand,  which  part  provides  more  contribution  to  calculation  of 
the  kd-pd  diagram. 

Since  in  (11),  i  takes  the  values  of  1  and  2,  we  have  a  set  of  four  homogenous  equations  to  be 
solved  for  deriving  the  desired  kd-pd  relation.  The  dispersion  relation  is  then  obtained  equating  the 
determinant  of  the  coefficients  to  zero  (note  that  the  variables  are  bto ,  b]^ ,  b'o  ^nd  ).  It  is  worth 
mentioning  that,  the  series  in  (14,  15)  are  slowly  convergent  series  and  hence  excessive  number  of 
terms  are  required  for  the  desired  preciseness.  Acceleration  techniques  have  been  applied  to 
numerically  solve  the  above  transcendental  equations  for  the  propagation  constant  p  [7,  8]. 

The  above  technique  is  successfully  applied  to  metamaterials  characterization.  Fig.  5(c)  obtains  the 
dispersion  diagram  for  two-set  of  spheres  having  the  same  size  and  different  materials.  The 
dispersion  diagram  and  the  existence  of  a  negative  slope  branch  are  demonstrated.  A  full  wave 
numerical  technique  based  on  Finite  Difference  Time  Domain  (FDTD)  technique  [9]  is  also  applied  to 
validate  the  accuracy  of  our  developed  method  based  on  dipole-mode  analysis.  An  excellent 
comparison  is  demonstrated,  where  our  method  is  much  more  efficient. 

An  array  of  spheres  with  the  unit-cells  made  of  two  different  spheres  (different  materials  and  same 
sizes)  is  modeled  in  Fig.  6.  The  design  is  optimized  such  that  the  electric  resonance  of  one  set  of 
spheres  will  be  around  the  magnetic  resonance  of  another  set.  Fig.  6(a)  and  6(b)  illustrate  the  electric 
and  magnetic  band-gaps  (for  each  set)  associated  with  the  relevant  resonances.  Combining  these  two 
sets  of  spheres  demonstrates  a  pass-band  through  the  stop  band  region  highlighting  a  backward 
metamaterial  structure.  One  can  also  design  a  negative-slope  backward  wave  medium  utilizing 
different  sizes  spheres  (having  the  same  materials)  in  a  unit  cell,  as  determined  in  Fig.  7.  It  must  be 
mentioned  that  the  backward  wave  is  occurred  in  the  region  that  kd  is  larger  than  1  and  effective 
materials  parameters  may  not  be  defined  successfully.  One  can  increase  the  dielectric  materials  of 
the  spheres  to  bring  the  backward  wave  dispersion  diagram  in  the  region  that  kd<l  and  thus  define 
properly  the  effective  material  parameters,  as  is  illustrated  in  Fig.  8.  This  realizes  a  metamaterial  with 
Double  Negative  (DNG)  materials  characteristics,  although  the  bandwidth  of  the  metamaterial  is 
relatively  small.  The  effect  of  the  coupling  on  the  dispersion  diagram  and  its  bandwidth  performance 
is  investigated  in  Fig.  9.  This  is  one  set  of  spheres  operating  in  their  magnetic  resonances  where  their 
couplings  provide  electric  modes.  The  hybrid  modes  tailor  a  dispersion  diagram  with  negative  slope. 
As  one  can  see  increasing  the  couplings  between  the  elements  accomplishes  more  bandwidth  for  the 
performance  which  will  be  a  very  novel  idea.  Full  investigations  and  discussions  about  the 
configurations  and  results  are  provided  in  [6]. 


Fig.  4:  Geometry  of  two-sets  of  dielectric  spheres  metamaterials. 


(a)  (b)  (c) 

Fig.  5:  Electric  and  magnetic  Mie  scattering  coefficients  of  a  single  sphere  for  the  first  three  modes  with  ,  (a) 
0  =  0.5cm  and  e^=AO,  (b)  o  =  0.5  cm  and  £*^=21,  (c)  Comparison  between  FDTD  results  and  theoretical 
approach  for  two-sets  of  dielectric  spheres  with  0^=03  =0.5 cm,  d/a  =  3,  hja  =  h^/a  =  5 ,  h  =  hJ2, 
and  =  21.  A  very  good  agreement  between  our  theoretical  model  and  FDTD  is  observed. 


Fig.  6:  Dispersion  diagrams  for  two-sets  of  dielectric  spheres  having  the  same  sizes  (Oj  =03  =0.5cm)  and 
different  dielectric  materials,  (a)  The  performance  of  the  first  set  described  with  ^^,i=40  and  c/  =  l.lcm, 
=2.5cm,  (b)  The  performance  of  the  second  set  described  with  6:^3  =20.5,  d  =  l.lcm  and 
=2.5  cm,  (c)  The  performance  of  the  combinations  of  the  two-sets  with  h^hjl . 


Fig.  7:  Dispersion  diagrams  for  two-sets  of  dielectric  spheres  having  the  same  dielectric  materials 
(^ri  =  ^r2  =  ^0)  different  radii,  =  0.5cm  and  =  0.356cm,  d  =  1.1cm  ,  =1.8 cm,  h-hj2  . 


Fig.  8:  Dispersion  diagrams  for  two-sets  of  dielectric  spheres  having  the  same  sizes  and  different  dielectric 
materials  =98,  =200,  Oj  =  0.5  cm,  d/a -2.1 ,  h^=hy  =2.5  cm  and  h-hj2  .  kd  <  1  is  achieved 

where  one  can  define  effective  materials  parameters. 


Fig.  9:  Dispersion  diagram  characteristic  for  the  second  branch  of  array  of  highly  coupled  cubical  unit-cell 
spheres  having  c^=20,  o  =  0.5cm  and  =/7^  =d  .  Increasing  the  coupling  will  change  the  performance  from 
positive  slope  to  negative  slope. 


3.  Magnetic  Materials-Based  Spheres  Metamaterials 


We  also  investigate  the  performance  of  metamaterials  realization  utilizing  array  of  magneto-dielectric 
spheres  or  dielectric  spheres  embedded  in  a  host  natural  hexa-ferrite  medium  operating  in  its 
negative  permeability  region  [10].  The  advantage  of  integrating  magnetic  materials  is  achieving  wider 
bandwidth  as  has  been  explored  earlier  by  Mosallaei  in  [11]. 

Our  developed  formulation  is  general,  and  Fig.  10  demonstrates  the  characteristics  for  an  array  of 
magneto-dielectric  spheres  [6].  The  scattering  coefficients  and  dispersion  diagram  are  shown  in  Fig. 
10.  The  scattering  coefficient  is  for  a  single  sphere  highlighting  the  electric  and  magnetic  resonances 
around  the  same  frequency  region.  Making  an  array  of  spheres  considering  the  couplings  between 
them  realizes  a  backward  wave  medium  in  the  region  that  kd<l,  and  where  a  wider  bandwidth  in 
compared  to  only  dielectric  design  is  established.  The  effective  materials  parameters  are  obtained  in 
Fig.  11  proving  a  DNG  medium  with  both  negative  s  and  p. 

We  also  exploit  a  novel  design  of  one  set  of  array  of  dielectric  spheres  embedded  in  a  natural  hexa- 
ferrite  material  with  Lorentzian  property  given  by  [10]  (and  plotted  in  Fig.  12(a)) 


(17) 


The  dielectric  particles  are  designed  to  operate  in  their  electric  resonances  (around  the  spectrum 
where  the  host  material  has  its  own  negative  permeability).  The  performance  is  determined  in  Fig. 
12(b).  The  effective  materials  associated  with  this  design  are  plotted  in  Fig.  13.  In  compared  to 
previous  design  a  better  performance  with  more  negative  permeability  range  is  developed. 


The  integration  of  magnetic  materials  can  also  allow  tunability  by  applying  magnetic  DC  fields. 
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Fig.  10:  (a)  MIe  scattering  coefficients  for  a  single  magneto  dielectric  sphere  with  =  5,  o  =  0,85cm  , 

(b)  dispersion  diagram  characteristic  for  one-set  of  magneto  dielectric  spheres  with  a  =  0,85  cm ,  d/a  =  2,1  and 
h. 


Fig.  11:  Effective  permittivity  and  permeability  for  3D  array  of  magneto-dielectric  spheres  with  cubical  unit  cell 
and  a  =  0.85cm  ,  =18,  =  5,  d/a  =  2.1 .  The  bulk  materials  are  plotted  In  the  region  where  both  kd^fdd  <  1, 


(a)  (b) 

Fig.  12:  (a)  Plot  of  the  Lorentzian  permeability  function  given  by  (17).  (b)  Dispersion  diagram  characteristic  for 
one-set  of  dielectric  spheres  embedded  in  a  Lorentzian  host  medium  with  cubical  unit-cell  and,  o  =  0.85cm, 
d/o  =  2.1,^^  =  90. 


0.987  0.988  0.989  0.99 

kd 


Fig.  13:  Effective  parameters  for  a  cubical  3D  array  of  dielectric  spheres  described  with,  s^=90,  o  =  0.85cm  , 
d/o  =  2.1,  inside  a  negative  permeability  material  with  Lorentzian  function  given  by  (17).  The  bulk  materials 
are  plotted  in  the  region  where  both  kd,j3d  <1. 


4.  FDTD  Analysis  and  Other-Shapes  Configurations 

Although  all  the  results  presented  so  far  are  for  the  metamaterials  realization  by  spheres  particles, 
other-shapes  particles  can  also  create  DNG  metamaterials,  structures  such  as  disks  or  rods,  achieving 
the  required  electric  and  magnetic  dipole  modes.  Here,  the  creation  of  a  metamaterial  made  from 
array  of  two  types  of  disks  is  investigated.  Characterization  is  based  on  a  full  wave  Finite  Difference 
Time  Domain  (FDTD)  technique. 

We  will  first  start  with  a  metamaterial  array  of  dielectric  spheres  and  then  we  will  move  to  the  array 
of  disks.  Fig.  14(a)  shows  the  geometry  and  effective  materials  parameters  for  the  array  of  spheres. 
The  FDTD  is  applied  and  transmission  coefficient  is  determined  in  Fig.  14(b).  A  band-pass  region  in  the 
DNG  band  is  achieved.  The  FDTD  transmission  coefficient  result  is  in  great  agreement  with  the 
dispersion  diagram  determined  by  our  dipole-mode  analysis.  The  near-fields  inside  the  spheres  are 
shown  in  Fig.  14(c).  It  is  clearly  observed  that  one  set  of  spheres  is  equivalent  to  an  electric  dipole  and 
the  other  set  equivalent  to  a  magnetic  dipole.  Combination  of  these  two  modes  features  a  backward 
wave  medium. 

Fig.  15  illustrates  a  metamaterial  made  of  array  of  disks  and  its  characteristics.  Electric  and  magnetic 
resonances  and  DNG  region  are  demonstrated.  The  near-fields  are  shown  in  Fig.  16,  clearly  presents 
the  electric  and  magnetic  modes  creation. 

A  full  discussion  on  FDTD  modeling  of  array  of  dielectric  particles  metamaterials  along  with  a 
comprehensive  numerical  results  and  physical  understandings  have  presented  in  [12]. 


Fig.  14:  FDTD  analysis  for  3D  array  of  spheres  providing  DNG  medium,  (a)  the  geometry,  (b)  transmission 
coefficient  and  dispersion  diagram,  and  (c)  near-field  patterns. 


(a)  (b) 

Fig.  15:  DNG  metamaterial  constructed  from  all-dielectric  disks,  (a)  the  geometry  ( =2.5cm,A^  =1.5cm), 
and  (b)  its  transmission  coefficient. 


Fig.  16:  Field  distributions  inside  one  unit-cell  of  the  DNG  metamaterial  (Fig.  15(a))  at  /  =  5.97 GHz  .  (a)  E  in  y-z 
plane,  and  (b)  H  in  y-x  plane.  Note  to  the  creation  of  electric  and  magnetic  dipole  moments  inside  the  unit-cell 
of  the  disks  of  8r=60  and  8r=43. 


5.  Magneto-Dielectric,  ENG,  and  MNG  Spheres  Small  Antennas 

This  section  explores  the  performance  characterization  of  small  antennas  embedded  in  magneto¬ 
dielectric,  Epsilon  Negative  (ENG),  and  Mu  Negative  (MNG)  spheres.  These  structures  can  miniaturize 
the  antennas  while  at  the  same  time  a  relatively  wide  impedance  bandwidth  can  be  established.  A 
full  detailed-study  can  be  found  in  [13-17].  Here,  we  summarize  the  findings. 


5.1  Meta  materia  Is  Antennas 

Fig.  17  depicts  the  geometry  of  an  antenna  embedded  in  a  hemi-sphere  metamaterial  located  on  a 
ground  plane.  The  resonator  of  radius  a  is  coupled  to  a  transmission  line  with  a  small  monopole  stub 
of  length  /  which  is  located  at  a  distance  b  from  the  center  of  sphere.  The  objective  is  to  obtain  the 
input  impedance  of  the  structure.  Using  the  equivalence  principle  and  image  theory  one  can  simply 
consider  the  structure  as  a  dipole  inside  a  sphere.  Obviously,  the  input  impedance  of  the  original 
problem  is  equal  to  one-half  of  the  new  structure. 


To  obtain  the  electric  field  in  the  z-direction  due  to  the  z-directed  current  element,  we  first 
decompose  the  current  vector  into  the  r-  and  0-components.  Note  that  r-directed  current  can  only 
generate  tm' modes,  where  the  0-directed  current  can  generate  both  TM'  and  TE'  modes.  For  the  0- 
directed  current,  we  represent  three  Green's  functions  as  incident  field  (dipole  only),  scattering  field 
inside  the  sphere  (the  field  inside  the  sphere  caused  by  the  sphere  discontinuity),  and  the  scattering 
field  outside  the  sphere;  as  given  below  [13]. 


"  ^  ^  [Hl^^ikr)  r>r 

(18a) 

n=0  m=-n 

(18b) 

n-0  m=-n 

(18c) 

Where,  j„(x),H^(x),P„'”(x)  are  spherical  harmonics.  The  unknown  coefficients  of  incident  and  scattered 
fields  are  obtained  by  applying  the  boundary  conditions.  e„,e^  and  Hg  are  continuous,  and  is 

discontinuous  by  0-directed  surface  current  distributed  on  r  =  r'  (note  that  ( r' ,0' )  refers  to  source 
point  while  {r,0,(f)  refers  to  field  point).  By  matching  the  boundary  condition  at  r=a  the  scattering 
solution  for  inside  and  outside  are  derived. 


The  next  step  is  to  derive  the  following  integral  equation  for  the  z-component  of  electric  field: 

f,  =  |j[(cos^Gj^'  -sin^Gj^')cos^  -(cos^G^;  -sin^Gj^')sin^']/^dS  (19) 

h 

Where  is  the  surface  on  which  the  current  is  flowing  and  Xis  the  z-directed  current.  Combining 
the  incident  field  and  the  scattering  field,  the  z-directed  Green's  function  due  to  the  z-directed 
current  can  be  written  as: 
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The  method  of  moment  with  Galerkin's  procedure  and  Piecewise  Sinusoidal  (PWS)  basis  functions  are 
applied  to  solve  for  the  current.  A  computer  program  is  developed  based  on  the  above  equations  to 
successfully  characterize  the  performance  of  an  antenna  embedded  inside  the  metamaterial  sphere, 
as  will  be  highlighted  next. 


The  performance  of  the  antenna  embedded  in  dielectric,  magneto-dielectric,  and  ENG  hemi-spheres 
are  shown  in  Figs.  18  and  19.  Both  ENG  and  magneto-dielectric  metamaterials  have  better 
performances  compared  to  only  dielectric  sphere  case.  It  is  shown  in  [11]  that  a  magneto-dielectric 
material  with  both  £  and  p  parameters  can  enhance  the  bandwidth  as  also  explored  here  (the 
matching  is  difficult  since  a  high  impedance  characteristic  is  obtained).  ENG  antenna  can  successfully 
miniaturize  the  antenna  size  [14-16].  From  equation  (21a),  considering  a  small  dimension  for  sphere, 
one  can  determine 
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Negative  permittivity  parameters  can  make  a  small  antenna  resonator  as  illustrated  in  Fig.  19.  The 
antenna  can  resonate  at  a  frequency  of  interest  almost  independent  of  the  size,  but  of  course  the  Q  is 
strongly  related  to  the  size  [18].  The  major  issue  is  that  one  might  not  be  able  to  naturally  obtain  a 
medium  with  negative  e  in  microwave  frequencies.  Magnetic  materials  are  more  promising  in  this 
regard  as  will  be  explored  in  next  section. 


5.2  MNG  Antennas 

One  can  use  a  self-biased  hexaferrite  in  GHz  spectrum  and  operates  above  the  material  resonance  to 
provide  desired  negative  permeability  [10].  Hence,  it  \will  be  of  practical  importance  if  one  can  design 
a  small  antenna  utilizing  MNG  materials  [17].  The  system  \will  be  dual  of  ENG  antenna,  and  one  needs 
to  excite  the  structure  \with  a  slot  excitation  realizing  magnetic  current.  Note  that,  the  slot  feeding 
system  is  more  convenient  for  fabrication  than  the  dipole  feed  that  must  be  drilled  into  the  ENG 
medium. 


Fig.  20(a)  depicts  a  hemispherical  resonator  of  radius  a,  composed  of  a  medium  \with  permittivity  of 
€  and  permeability  of  and  excited  by  a  slot  aperture  of  length  21  and  width  of  W  (providing  y- 
directed  magnetic  dipole).  To  successfully  present  a  Green's  function  analysis  of  the  structure,  the 
electromagnetic  fields  are  decomposed  into  the  contributions  from  the  stripline  and  the  slot 
excitation,  and  the  scattered  field  due  to  the  sphere  discontinuity  [17].  Applying  the  reciprocity 
theorem  along  with  the  assumption  that  the  stripline  has  an  infinite  length  allow  representing  the 
antenna  configuration  by  the  following  equation: 

K]  =  {[C  +C]+^[A'/J[AvJ j  [AvJ,  (23) 


where  contributes  for  the  effect  of  the  slot,  is  the  admittance  caused  by  the  discontinuity  of 
the  sphere,  and  Av  is  associated  with  the  slot's  discontinuity  voltage  due  to  the  stripline  excitation. 
Y^  is  obtained  by  deriving  the  Green's  function  which  is  the  Fourier  transform  of  at  (x,y,0) 
due  to  a  unit  y  -magnetic  current  located  at  (0,0,0).  It  is  given  by  [19]: 
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One  can  determine  Av  by  deriving  the  Green's  function  6™,  which  is  the  Fourier  transform  of  at 
(x,y,0)  due  to  a  unit  x-electric  current  element  at  (0,0,z/),  it  is  specified  by: 
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The  admittance  Y^  is  resulted  from  the  sphere  discontinuity  and  has  direct  impact  on  the  resonant 
frequencies.  To  evaluate  it,  the  Green's  function  for  a  y-directed  magnetic  current  located  at  the 
center  of  the  sphere  should  be  obtained.  It  is  determined  that  [17] 
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where: 


b„=(2n  +  l)A 


(27) 


The  characteristic  equation  is  then: 
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For  characterizing  small  size  sphere  with  ka.k^a  «0.5  it  would  be  very  beneficial  if  one  can  obtain  an 
approximated-simple  formulation  for  the  Green's  function  along  the  slot  (y-axis).  For  this  case,  one 
can  simplify  (approximate)  equation  (26)  as  follow: 
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The  resonant  frequency  of  MNG  antenna  can  be  predicted  from  the  characteristic  equation 
when  =  0,  where 

^l^j{klk,y  ^  (f?//^+f?  +  l),  (31) 

(2n  +  l)// 

It  is  observed  that  the  antenna  can  resonate  for  n,  =-(n  +  l)/n ,  which  is  interestingly  independent  of 
the  sphere  size.  The  first  dominant  mode  (n  =  l)  offers  resonant  frequency  at  n,=-1  (quasi-static 
model).  Higher  order  modes  resonate  at  other  negative  permeability  values. 


Using  (24),  (25)  and  (26)  (or  (29))  one  can  derive  and  Y^^  in  an  integral  equation  form  depending 
on  the  MoM  basis  functions  of  the  voltage  across  the  slot  [20].  The  matrix  equation  (23)  can  be 
solved  to  obtain  [\/J,  and  from  that  calculate  series  slot  impedance  Z,  as: 
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where  is  the  characteristic  impedance  of  the  stripline.  The  is  tuned  by  the  open  circuit  stub 


length  L^,  providing  total  antenna  input  impedance: 

Z,„=Z,-jcot(j3,L^),  (33) 

and  Pf  is  the  propagation  constant  in  the  stripline. 


To  validate  the  applicability  of  the  developed  formulations,  the  performance  of  a  MNG  hemisphere 
constructed  from  a  Lorentzian  medium  is  investigated.  The  hemisphere  has  radius  of  o  =  7.5mm  and 
is  excited  with  slot  feed  having  length  of  2/  =  12mm  and  width  of  W  =  0.9mm.  The  Lorentzian 
permeability  function  is  given  by  (plotted  in  Fig.  20(b)) 

2  6)^  -  JIttxAxIO^ 6) ^ 


1-- 


(34) 


Numerical  results  for  the  input  impedance  based  on  approximated  and  exact  solutions  (equations 
(29)  and  (26))  are  presented  in  Fig.  21,  verifying  the  accuracy  of  the  developed  formulations.  The 
resonant  frequencies  are  around  f  =  2.18GHz  and  f  =  2.27GHz,  associated  with  ^^=-2.09 
(fundamental  mode)  and  //,  =-1.23  (third  mode),  respectively.  These  results  successfully  validate  the 
applicability  of  equation  (31)  in  predicting  the  resonant  frequencies.  As  observed  in  Fig.  21,  the  first 
dominant  mode  can  provide  a  wider  impedance  bandwidth.  The  return  loss  is  shown  in  Fig.  22,  where 
it  demonstrates  antenna  impedance  matching  at  the  first  resonant  frequency.  The  -lOdB  bandwidth 
is  about  0.45% .  Tuning  the  stub  length  properly  matches  the  antenna  at  the  first  resonant  frequency. 
The  total  size  of  sphere  at  the  fundamental  mode  is  smaller  than  >1/9  ( ka  =  0.34 ). 


Fig.  17:  Configuration  of  a  metamaterial-based  antenna. 


Fig.  18:  Input  impedance  and  return  loss  vs.  frequency  using  delta  gap  source  model,  (a)  Input  impedance,  (b) 
Return  loss:  a  =  12.5  mm,  b  =  6.4  mm,  /  =  6.5  mm,  =  9.8,  =  0.063  mm. 


Re*l  Part  Imaginary  Part 


Fig.  19:  Input  impedance  for  ENG  and  magneto-dielectric  antennas,  (a)  Input  resistance  (b)  Input  reactance:  a  = 
7.5  mm,  b  =  0,  /  =  6  mm,  0.4  mm.  Drude  medium  with  4  =46Hz,^^  =0.00lrj, ,  and  magneto-dielectric 
material  with  =7.8,£r,  =6.73). 
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Fig.  20:  (a)  The  geometry  of  MNG  hemisphere  antenna  with  radius  of  a  =  l,5mm  located  above  an  Infinite 
ground  plane,  fed  by  slot  with  length  of  2/  =  \2mm  and  width  of  W=0.9  mm.  The  stripline  has  dimensions  of 
Wj.  =\A5mm  and  L,  ^  2.6cm,  The  substrate  has  permittivity  of  =2.96  and  thickness  of  d  -0.653mm .  (b) 
The  plot  of  Lorentzlan  permeability  function  given  by  (34). 


Fig.  21:  Input  impedance  of  the  MNG  hemisphere  antenna  obtained  using  the  approximated  and  exact  Green's 
functions  analysis,  (a)  Input  resistance,  and  (b)  input  reactance. 


Fig.  22:  Return  loss  of  the  MNG  hemisphere  antenna. 
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